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We measured the production of the long-Hved nucHdes ^"^Bi, ^"^Pb, and '^^Hg in a sample of natural Pb due to high-energy neutron 
interactions using a neutron beam at the Los Alamos Neutron Science Center The activated sample was counted by a HPGe detector 
to measure the amount of radioactive nuclides present. These nuclides are critical in understanding potential backgrounds in low 
background experiments utilizing large amounts of Pb shielding due to cosmogenic neutron interactions in the Pb while residing on 
the Earth's surface. By scaling the LANSCE neutron flux to a cosmic neutron flux, we measure the sea level cosmic ray production 
rates of 8.0 + 1.3 atoms/kg/day of ''*'*Hg, 120 ± 25 atoms/kg/day ^"^Pb, and 0.17 + 0.04 atoms/kg/day ^''^Bi. 
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■1. Introduction 

I Neutrinoless double-beta decay iOv/3/3) plays a key role in 
understanding the neutrino's absolute mass scale and particle- 
antiparticle nature fHQ-SBSSBl- If this nuclear decay 
process exists, one would observe a mono-energetic line origi- 
nating from a material containing an isotope subject to this de- 
cay mode. The key to these experiments lies in the ability to 
reduce intrinsic radioactive background to unprecedented lev- 
lels and to adequately shield the detectors from external sources 
of radioactivity. Previous experiments' limiting backgrounds 
have been trace levels of natural decay chain isotopes within 
the detector and shielding components. The y-ray emissions 
from these isotopes can deposit energy in the detectors produc- 
ing a continuum, which may overwhelm the potential Ovp/3 sig- 
.nal peak. Great progress has been made in identifying the lo- 
'cation and origin of this contamination, and future efforts will 
substantially reduce this contribution to the background. The 
background level goal of 1 event/ton-year, however, is an ambi- 
tious improvement over the currently best achieved background 
level 1 8, 9, 10]. If the efforts to reduce the natural decay chain 
isotopes are successful, previously unimportant components of 
the background must be understood and eliminated. The con- 
tribution from long-lived isotopes produced by cosmic-ray neu- 
trons in detector and shielding materials must be considered. 
For example, the work of Ref. fll'] measured cosmic activa- 
tion of ^^Ge. Several Qv/3/3 experiments utilize large amounts of 
low background Pb as a means to passively shield external ra- 
diation. If the Pb shielding resides on the Earth's surface prior 
to deployment at an underground laboratory, cosmic activation 
occurs. No measurements of cosmic activation in Pb exists or 
lines identified as such by low background experiments. It is 
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not clear if the lack of prior observations is due to short surface 
exposures of low background Pb or due to a small activation 
rate. Therefore, it is important to understand the production of 
long-lived activated products in Pb for future Ov/?/? decay exper- 
iments 

We exposed a sample of Pb to a wide-band neutron beam 
that resembles the cosmic -ray neutron flux. After exposure we 
counted the sample in a low -background counting system to ob- 
serve the 7 rays from the decays of problematic isotopes. From 
these data we measure the production rate due to fast neutrons 
in the Pb sample. With knowledge of the neutron-beam and 
cosmic-neutron energy spectra, we use these data to provide an 
estimate of the production rate due to exposure of Pb to cosmic 
rays. We use a cross-section calculation that spans the energy 
range of interest to validate a comparison of production due to 
a neutron beam and the most recent cosmic -ray neutron flux 
measurements of which we are aware. This article describes 
our determination of the production rate of these isotopes. 

2. Experiment 

A natural Pb sample was exposed to the neutron beam at 
the Los Alamos Neutron Science Center (LANSCE) Weapons 
Neutron Research (WNR) facility from Target 4 Flight Path 
60 Right (4FP60R) fl^. As the broad-spectrum, pulsed neu- 
tron beam strikes the Pb target, the outgoing y rays are de- 
tected by the GErmanium Array for Neutron Induced Excita- 
tions (GEANIE) spectrometer 01 311 . The corresponding data 
from the GEANIE spectrometer and (n,n'7) analysis was pre- 
sented in a separate publication lfl4ll . The GEANIE sample is 
located a distance of 20.34 m from the natural tungsten spalla- 
tion target. 

The neutron target at the center of GEANIE was five stacked 
foils of natural Pb ('^^'■Pb) angled 20° from the normal to the 
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Figure 1 : The energy spectrum of the neutron beam at 4FP60R for the 2006 
exposure periods. The solid curve represents the sea level cosmic-neutron flux 
from Ref. fUl. 



Figure 2: The production cross sections for ^"^Pb (top curve) and "^Hg (lower 
curve) calculated by TALYS for neutrons impinging on natural Pb. The neutron 
energy range of TALYS is limited to 250 MeV. 



beam direction. Each foil measured nominally 5 cm x 5 cm 
in area and 0.475 mm in thickness. The neutron exposure oc- 
curred over two periods with 2 foils exposed in the year 2003 
and all 5 foils later in year 2006. The first 2-foil exposure was 
performed between Aug. 18 and Aug. 22, 2003 (3.72 days). 
The full 5-foil exposure occurred between Oct. 13 and Oct. 30, 
2006 (16.70 d), Nov. 9 and Nov. 13, 2006 (3.66 d), and on Nov. 
22 (0.33 d). The pulsed neutron beam has the following tim- 
ing structure. Sub-nanosecond-long neutron micropulses occur 
every 1.8 fis during a 625 //s-long macropulse, which occurs at 
a nominal rate of 40 Hz. The neutron energy is determined by 
the time of flight from the micropulse start. An in-beam fission 
chamber measures the neutron flux with ^■'^U foils |15]. As 
seen in Fig.[Tl the neutron energy spectrum at 4FP60R in 2006 
shows good spectral agreement over an energy range between 
20-300 MeV with the Gordon 1 16] parameterization of the sea 
level cosmic-neutron energy spectra. If the reader wishes to 
convolve the neutron spectrum with his/her own cross section 
model, we give a parameterization of the neutron spectrum im- 
pinging upon our sample for convenience. The spectrum can be 
described as: 



<!)(£) = (1.832 X 10'°)x (1) 

g(5.1661n£-3.815 In^ £+0.895 In' £-0.071 In* £) 

where O is in units of neutrons/MeV and the energy (E) is in 
MeV. 

The TALYS UtIi nuclear reaction code was run to to calcu- 
late excitation functions for producing nuclides in Pb. Of the 
nuclides predicted by TALYS, only '^"^Hg and ^"^Pb are long- 
lived (fi/2 > 1 yr) and support radioisotopes with a decay mode 
to an excited state so that y-ray emission is possible. '^"^Hg sup- 
ports the shorter lived ^^"^Au (fi/2 = 38.0 h) and ^''^Pb supports 
the shorter lived 202^] ^^^^^ - 12.3 d) 11811 . The production cross 
sections for these two nuclides are shown in Fig. |2l The pro- 
duction cross section of ''''*Hg includes feeding from the short- 



lived nuclides '^"^Tl and '^'*Pb. Note that the cross sections are 
peaked within the energy range of 20-300 MeV, which is the 
same energy region where there is agreement in the shape of the 
neutron flux between sea level cosmic neutrons and the 4FP60R 
neutrons (Fig. [T]). 

The experimental setup at GEANIE is optimized for mea- 
suring prompt y rays due to neutron reactions in a target sam- 
ple. Hence, the neutron irradiation during periods while the 
data acquisition system is not running does not affect prompt 
in,n'y) analysis but could constitute an appreciable neutron flu- 
ence for an activation study. Upstream of the WNR target exists 
a proton current monitor that continuously logs the beam cur- 
rent directed at WNR. This monitor provides a better handle 
on the total neutron current delivered to the flight path if the 
fission chamber neutron monitor experiences periods of down- 
time (e.g. in between data runs, during DAQ malfunctions and 
debugging, etc.). A subset of the irradiation period while the 
neutron fission chamber is online can be used to normalize the 
proton current to delivered neutron current. 

After irradiation, the Pb foils were stored for an extended 
period and therefore any radioactivity had decayed to an ex- 
tremely low level (<10 Bq) before counting began. Therefore, 
this sample was well below any action levels and not subject 
to any source-handling requirements. The sample was trans- 
ported to our low-background counting facility underground at 
the Waste Isolation Pilot Plant (WIPP) near Carlsbad, NM and 
counted on a HPGe detector. The detector was fabricated in 
1985 and placed underground at WIPP in 1998. It is an n-type 
semi-coax design with a height of 41 mm and a diameter of 
51 mm. It is contained within an «il-mm thick Cu cryostat. 
The shield during these runs consisted of 5 cm of oxygen-free, 
high-conductivity Cu and 10 cm of Pb. The detector and its Pb 
shielding has been underground since 1998. 

The sample was counted over a period of 126.6 days between 
Aug. 4 and Dec. 8, 2010 with a total live time of 109.07 days. 
The y ray energy spectrum is shown in Fig. [3] Three lines from 



2 




100 200 300 400 500 600 700 800 

Energy [keV] 



In addition to the y-ray detection efficiency (e^), an efficiency 
factor must also be included to take into account the decay of 
the isotope since the end of exposure and during the counting 
period. This latter efficiency factor (e^) depends on the half-life 
of the isotope of interest and is also given in Table [1] 

With these efficiencies, it is straight-forward to calculate the 
number of atoms of each isotope that were present at the end 
of neutron exposure (November 22, 2006) and the production 
rate. 

The measured number of counts (C) is related to the number 
of atoms (A^,) on the reference date at the end of exposure by: 



C 



(2) 
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Figure 3: The energy spectrum of y rays from the Pb sample (red) as mea- 
sured by a Ge detector and a background spectrum (blue) taken with no sample 
present. 



"^Au, a line from ^"^Tl, and two lines from ^"^Bi are observed 
in the activated spectra in addition to background lines. 

3. Measured Production Rates 

The peaks in Fig. |3] were fit to determine the measured counts 
(C) and the results are given in Table [T] In some of the peaks, 
the same lines are observed in the background spectra. In the 
case of ^°^Bi, the 569-keV and 1064-keV y ray lines are also 
found in the background spectra at a lower rate. The 295-keV 
y-ray line attributed to ''^"^Au is also found in the background 
spectra due to naturally occurring ^''^Pb. The measured counts 
in Table[T]represent the net counts after background subtraction 
for these three lines. Note that the 328-keV line could be due to 
decay of naturally occurring ^^'^Ac. However, the background 
spectrum shows no evidence of this Ac line and therefore the 
entire line is attributed to the decay of '^"^Au. 

The counting efficiency was determined by Monte Carlo 
simulations using the GEANT4-based MaGe framework 1 19.1. 
Simulations of the detector's counting efficiency is bench- 
marked against measurements of calibrated radioactive stan- 
dards (^^'^''Co, ^"^Mn, ^^Na and '■'^Cs). We use a random pulser 
to verify the event-rate dependence of the data-acquisition sys- 
tem dead time. 



where the sum is over the four irradiation periods. Since the 
first irradiation period only included two of the Pb foils, atoms 
produced from that run have a different y-ray detection effi- 
ciency than those produced in the latter three runs with all five 
Pb foils. The decay of the isotope before and during counting 
is corrected by 



-ATf 



-at; 



(3) 



where TJ""' 



(Tj""'') is the number of days since the end of ex- 



posure that the counting stopped (started) for each of the j data 
runs, and A is the decay constant of the isotope in question. 

The number of atoms can also be expressed in terms of the 
reaction rates during an exposure period: 



Ni 



(4) 



where we have corrected for the decay during the exposure and 
the decay after the exposure until the reference date. In Eqn.|4l 



decay 



jirraa duration of irradiation of the i exposure, and T\ 
is the time between the end of exposure / and the reference date. 

is the production rate (atoms/kg/day) during the LANSCE 
exposure /, and M, is the mass of the sample within the beam 
spot. 

The production rate of the identified long-lived nuclides in Pb 
occurs in an energy range where the shapes of the LANSCE and 
cosmic neutron fluxes are similar. Hence, a simple scaling of 
the LANSCE production rate is sufficient to estimate the cosmic 
production rate, given by 



K. 



irrad 



(5) 
(6) 



where K'^ is the cosmic production rate (atoms/kg/day), 
(/^) is the energy-integrated neutron flux for cosmic (LAN- 
SCE) neutrons over the mutual energy range 20-300 MeV, S 
is the normalization between proton current and neutron cur- 
rent, and A, is the cross-sectional area of the Pb target in the 
neutron beam. The sea-level cosmic neutron flux of Ref. lfl6|] 
is assumed where - 228 n/cm^/day. 
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Table 1: A summary of the measured cosmic production rates {K'') and the predicted cosmic production rates (^^^^lys-* ^™ long-lived isotopes in the Pb sample. 
Due to multiple y lines, we report the emphasized rate values given in column 8. If the observed nuclide is short lived, the half-Hfe of the supporting parent (listed 
in parentheses) is provided. The stated uncertainly on the production rate includes all statistical and systematic uncertainties combined in quadrature. 



Isotope 


T 1/2 


7 ray 


C 






ec 






"^TALYS 




[days] 


[keV] 




2006 


2003 






[/kg/day] 


[/kg/day] 




1.90x10^ 


293.5 


105+114 


0.0041+0.0002 


0.0042+0.0015 


3.95x10- 


-4 


7.0+7.7 






1.90x10^ 


328.5 


842+52 


0.0287+ 0.0014 


0.0294+0.0105 


3.95x10- 


-4 


8.0+1.3 


16 




1.90x10^ 


621.7 


101+18 


0.0009+0.00005 


0.0010+0.0003 


3.95x10- 


-4 


29+7 




202^1 (202pb) 


1.93x10^ 


439.5 


188+28 


0.0444+ 0.0022 


0.0453+0.0108 


3.90x10 


-6 


120+25 


77 


2'"Bi 


l.lTxlO'^ 


567.7 


316+37 


0.0341+0.0017 


0.0346+0.0055 


5.94x10 


-3 


0.17+0.03 




207Bj 


1.17x104 


1063.7 


146+26 


0.0153+ 0.0008 


0.0153+0.0011 


5.94x10- 


-3 


0.17+0.04 





Eqn.|4]can be written as 

/ irrad 

where the Mi/Ai term is replaced by the areal density (n,) of the 
Pb target parallel to the beam direction. The beam spot covers 
most, but not all, of the sample geometry. 

Table [T]lists the isotopes found and the number of measured 
counts after irradiations. Eqn. |7] allows a direct determination 
of the cosmic production rate K*-, also listed in Table [T] Of 
the three '^"^Au lines observed, two have low branching ratios 
and so we only place the most confidence on the 238-keV line. 
Further, the 622-keV line can occur as a lone y ray or as a 
sum of the 293.5-keV and 328.5-keV y rays. Its excess may 
be due to an additional systematic in Monte Carlo of the y ray 
efficiency. Further, there may be presence of the 622-keV y 
ray from '"*Rh, which is progeny of fission product '"*Ru pre- 
dicted by TALYS to be activated at a rate of 0.05 atoms/kg/day. 
If we assume the excess of the 622-keV line is due to '*'*Rh, 
we find a cosmic activation rate of 0. 18+0.06 atoms/kg/day for 
that isotope. The two ^"^Bi lines are found to have consistent 
production rates. 

Also hsted is the cosmic activation rate calculated with the 
TALYS excitation functions (Fig. |2]and the sea level cosmic 
neutrons [16]). The production of ^*"Bi in Pb requires an in- 
crease in the atomic number from secondary reactions from 
incident neutrons. Hence, the production is not predicted by 
TALYS. 

The aerial densities and total neutron yield for the two irra- 
diation periods is listed in Table |2] The low neutron yield and 
lower Pb thickness in the 2003 run make it a small contribution 
to the total production of nuclides. 

Table 2: A summary of the total neutron yield for the two iiTadiation years. 



Year 


n 


Neutrons on target (S x I) 




(g/cm^) 


between 20-300 MeV 


2003 


1.215 


5.73+5.46 X 10'" 


2006 


3.038 


6.77+0.28 X 10" 



The uncertainty in the production rate is dominated by the 
counting statistics of weak peaks (Table [T]i, the counting effi- 
ciency, the exposed neutron yield, and the cosmic ray neutron 



flux. The y-ray efficiency for the full five Pb foils exposed in 
2006 is known to 5%, which is based on source data and Monte 
Carlo agreement. The placement of the two Pb foils exposed in 
2003 within the five foils during y counting is unknown, though 
they were adjacent. Therefore, the 2003 y-ray efficiency has 
a larger systematic uncertainty (Table [TJ. The uncertainty on 
neutron yield is due to the statistical uncertainty on the proton 
beam current and the uncertainty on the proton to neutron nor- 
malization. The 2006 total neutron yield is known to 4.2%. For 
the 2003 neutron yield calculation, the limited fission chamber 
data available suggests a lower proton to neutron normaliza- 
tion, though no physical explanation can explain the discrep- 
ancy with the 2006 normalization. We assume the nominal 
2006 normalization for the 2003 neutron yield, but conserva- 
tively include a larger uncertainty of 95% to remain consistent 
with the limited data. The 2003 exposure represents 3.3% of 
the total neutron exposure so the large uncertainty on the 2003 
neutron yield has a small effect on the net production rate. 

The precision to which the cosmic ray neutron flux is known 
is 10-15% and so we split the difference and use 12.5% ifieli . 
The total cosmogenic rate includes contributions from sub- 
dominant proton and pion interactions. These only contribute 
approximately 10% ll20ll to the total production rate. These 
charged particles are much less penetrating than neutrons and 
therefore their impact on any given sample is very geometry 
dependent. Hence we assume a 50% uncertainly on this cor- 
rection for an uncertainty of 5%. Again these uncertainties are 
uncorrected and result in a estimated systematic uncertainty on 
the cosmic neutron flux of 13.5%. The start and stop times of 
counting and the live time of the counting are known to a small 
percentage and are negligible contributions to the uncertainty. 
This is similar for the times associated with the irradiation. The 
sample had been stored on the Earth's surface for many years 
prior to exposure to the beam and then counting at WIPP. Any 
isotopes produced by cosmic ray neutrons would certainly be 
below saturation after this extended period with a cosmic con- 
tribution •« 3% . The saturation production rates however are 
expected to be low. Therefore the total count rate is expected to 
be dominated by the intense neutron beam and we ignore this 
systematic effect. The values for the half-lives and branching 
ratios are known to high precision and are a negligible contri- 
bution to the total uncertainty. 
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4. Discussion and Conclusion 

We measured the production of '^"^Hg, ^°^Pb, and ^*"Bi in a 
sample of Pb due to high-energy neutron interactions within a 
neutron beam with a spectrum similar to that of the cosmic -ray 
neutron flux at the Earth's surface. These results can be used 
to predict sea level neutron activation of long-lived nuclides in 
Pb shielding. Previous low background experiments using Pb 
shielding have not found lines in their detectors energy spec- 
trum attributed to long-lived activated products in Pb. For ex- 
ample, the HEILDELBERG-MOSCOW (H-M) ^ neutrino- 
less double-beta decays experiment does not present evidence 
of '^'^Au or ^"^Tl. The H-M spectra does show the presence of 
^''^Bi, though its attributed to anthropogenic rather that cosmic 
production in Pb shielding. The ^"^Bi may be partially cosmo- 
genic in origin in the Pb shielding. If cosmogenic, the decay 
rate of ^"^Bi observed here suggests that ''"^Au and 202^] would 
also be observed. However, the activation of ^"^Bi here is due to 
secondary reactions involving protons and the net cosmic pro- 
duction would be greater due to the flux of cosmic protons. 

The work of Mei et al. fl^ using a Pb shield surrounding a 
HPGe detector did observe a line at 438.9 keV, which is suspi- 
ciously close the the 202'pj jjjjg 439.5 keV observed here. A 
long term spectra obtained by Ref. (2^ also observed a line at 
438.8 keV but attributed it to a "^''K double escape peak. The two 
also see the 328 keV Une from ^^^Ac, which could be blending 
with cosmogenic '^'^Au. However, the rate of the 328 keV line 
is consistent with the nearby ^^^Ac 338 keV line when scaled 
by their branching ratios. 

Neutrinoless double-beta decay experiments typically 
choose target nuclei with a high (> 2 MeV) double-beta Q 
value in order to stay above most y rays from natural sources. 
The activated products discovered here do pose a potential 
background contribution due to their high decay Q values 
and high energy y rays. While '^"^Au has a decay Q-value of 
2501 keV with y rays up to 2413 keV and ^'^^Bi has a decay 
Q-value of 2398 keV with y rays up to 1770 keV, ^"^ti is less 
of a concern with a lower Q value of 1363 keV lll8ll . The 



production rates of these cosmically activated radioisotopes in 
Pb can be used to estimate background contributions and sur- 
face exposure limitations for the next generation neutrinoless 
double-beta decay experiments. 

With its high energy y rays, '^"^Au is likely the most problem- 
atic cosmogenic in Pb shielding. The highest intensity y-ray 
line above 2 MeV occurs at 2043.7 keV. This transition is close 
to the 2039-keV Q value of double-beta decay target ^^Ge. Us- 
ing the production rate measured here, cosmogenic ''''*Au will 
saturate at a concentration of 93 /iBq/kg. With a decay intensity 
of 3.81%, the 2044-keV y ray would be produced at a rate of 
3.5 /^Bq/kg at saturation. As stated earlier, the H-M experiment 
did not see evidence of ''^'*Au. To estimate their sensitivity to a 
2 MeV y ray, they did observe the 2.6-MeV y ray from ^os-pj jjj 
the ^^^Th decay chain and determined the concentration to be 
12.3 yuBq/kg of ^^^Th (4 yuBq/kg of ^'^xi) in the Pb shielding 
1 2411 . Therefore, they would have sensitivity to ''^'*Au at satu- 
ration indicating the Pb shielding had little exposure to cosmic 
neutrons at the surface of the Earth prior to use underground. 



Current and future experiments are finding shielding materials 
with improved contamination levels. Pb suitable for shielding 
has been obtained at <4 y^Bq/kg of^^^Th (<1.3 A^Bq/kg of ^'"^Tl) 
1 2511 . Since ^^^Th contamination in Pb is currently unavoidable, 
it can be used to set the scale of other contaminate goals. In 
order to keep the decay rate of the 2044-keV y ray at or below 
the decay rate of the ^o^jj 2614-keV y ray, Pb surface exposure 
should be Umited to 350 yr A limitation of 30 yr of surface 
exposure is required to stay a factor of 10 less than the ^('S'pj 
decay rate. 
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